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Abstract: We analyze the phenomenology of a set of minimal supersymmetric standard 
model (MSSM) soft terms inspired by flux-induced supersymmetry (SUSY)-breaking in 
Type IIB string orientifolds. The scheme is extremely constrained with essentially only 
two free mass parameters: a parameter M, which sets the scale of soft terms, and the [i 
parameter. After imposing consistent radiative electro-weak symmetry breaking (EWSB) 
the model depends upon one mass parameter (say, M). In spite of being so constrained 
one finds consistency with EWSB conditions. We demonstrate that those conditions have 
two solutions for \i < 0, and none for fi > 0. The parameter tan/3 results as a prediction 
and is approximately 3-5 for one solution, and 25-40 for the other, depending upon M 
and the top mass. We examine further constraints on the model coming from b — * sj, 
the muon g — 2, Higgs mass limits and WMAP constraints on dark matter. The MSSM 
spectrum is predicted in terms of the single free parameter M. The low tan/3 branch is 
consistent with a relatively light spectrum although it is compatible with standard cos- 
mology only if the lightest neutralino is unstable. The high tan (5 branch is compatible 
with all phenomenological constraints, but has quite a heavy spectrum. We argue that 
the fine-tuning associated to this heavy spectrum would be substantially ameliorated if an 
additional relationship ^ = —2M were present in the underlying theory. 
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1. Introduction 

The MSSM is one of the most promising candidates for an extension of the Standard 
Model (SM). In applications to phenomenology a crucial ingredient is that of the structure 
of SUSY-breaking soft terms. A large number of models and scenarios for soft terms have 
been proposed in the literature (see e.g. ref.JTJ for a recent review and references). Some 
of the most popular schemes, like dilaton/modulus dominance SUSY-breaking and their 
generalisations are inspired by heterotic string models (see e.g.[|| for details and references). 

In this context it has been realised in the last few years that strings other than the 
heterotic, particularly Type II and Type I string theories, are equally viable as candidates 
in which to embed the MSSM. A first analysis of the soft terms in this class of models was 
presented in ref.||]. More recently it has been found that fluxes of antisymmetric 3-forms 
which are present in Type IIB string theory are natural sources for SUSY-breaking in Type 
IIB orientifold models 0|. In the latter the SM fields are assumed to live on the world- 
volume of D7 or/and D3 branes and/or their intersections. Specifically, if one assumes 
that the SM fields correspond to 'geometric moduli' of D7-branes, it was pointed out in 
ref.[|| that a very simple set of SUSY-breaking soft terms appear if certain background 
fluxes are present 1 . In fact these soft terms may be understood as coming from a modulus- 
dominance scheme applied to the particular case of Type IIB orientifolds, which leads to 
results different to those in the heterotic case. Although no specific string model with these 
characteristics has been constructed, the structure of soft terms is so simple and predictive 
that it is certainly worthwhile examining its phenomenological viability. 

For other structures of MSSM soft terms corresponding to a different localisation of SM particles on 
the D-branes see |], f| |. 
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The boundary conditions for the SUSY breaking terms are a subset of those of the 
minimal supergravity (mSUGRA) model. All scalars have a common mass mo, the gauginos 
have a common mass Mx/2 and the trilinear scalar couplings are identical to each other 
(once divided by the corresponding Yukawa coupling) and denoted Aq. In the mentioned 
scheme || such parameters are constrained as follows: 

M 1/2 =M, A = -3M, m = \M\ (1.1) 

where M parametrises the overall SUSY breaking mass scale in the model. In the simplest 
case M coincides with the gravitino mass. This would imply that the gravitino is not 
the lightest supersymmetric particle, so it tends to decouple from phenomenology, being 
very weakly coupled to matter. We should bear in mind, however, that the relationship 
between M and the gravitino mass may vary in particular models. Another independent 
mass parameter of the model is fj,, which appears in the superpotential term — iaH\H2. We 
parametrise the associated soft breaking term in the scalar potential as —iiBH\H%. Also 
the B soft parameter is predicted in terms of M: 

B = -2M (1.2) 



One of the nice features of the set of soft terms in eqs. (IT), (|l.2| ) is that complex phases 
may be rotated away, hence there is no 'SUSY-CP problem'. The MSSM with the above 
structure of soft terms was named as the 'fluxed MSSM' in || (where \iB was denoted as 
-B). 

Note that the above set of soft terms is extremely restrictive. There are only two free 
mass parameters fi and M. Imposing appropriate electroweak symmetry breaking (EWSB) 
essentially will lead us to a single parameter (say, M) corresponding to the overall scale 
of SUSY-breaking. Thus it is not at all obvious that consistent radiative EWSB may be 
obtained in such a constrained system. Remarkably, we find that the above soft terms are 
consistent with EWSB with tan/3 being fixed around two regions, with tan/3 = 3 — 5 and 
tan (5 = 25 — 40 respectively. The complete sparticle spectrum is then fixed depending only 
on M for those two regions. In what follows we will carry out a detailed analysis of the 
conditions of radiative EWSB in this scheme. We will also present constraints coming from 
b — > S7, the muon g — 2, Higgs mass limits and dark matter relic density. 

It is interesting to compare these results with those coming from other string-inspired 
schemes with a reduced number of soft parameters. One of the most attractive schemes is 
that of the dilaton-domination scenario |]] in heterotic string models in which the boundary 
conditions are given by M\/2 = —Aq = ±v3mo, mo being a universal scalar mass. In this 
case there are three free parameters corresponding to mo, /i and B, hence the general dilaton 
domination scenario is less predictive than the fluxed MSSM model here analysed which 
has only two free parameters, M and jx, B being fixed by eq. ( |1.2j ). It is thus particularly 
remarkable that correct radiative EWSB may be achieved in this very constrained system. 
There are, however, some restricted versions of the dilaton domination scenario which are 
equally constrained, and proper EWSB can be obtained even there. This happens, for 
instance, if jx is generated through a Giudice-Masiero mechanism [10], which leads to the 
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the further constraint B = 2itiq || 11]. Proper EWSB can be obtained 2 provided B and 



M1/2 have opposite sign [13], as in our case (see eq. Ql.2j) ). 

A few comments are in order about our notation and procedure. In general, we will 
write the soft terms in the notation of SDFTSUSY [14]] (except for fiB = m 2 ). We can take 
M real and positive without loss of generality. In order to discuss EWSB, it is convenient 
to replace the boundary condition eq. (|1.2[) with a more general form, that is 



B = -2Mr B . (1.3) 



The default case corresponds to re = 1. Eqs. (b_l), ( |Q| ) are to be applied at some high 



fundamental scale. For simplicity, we take such a scale to be Mqut, which is defined by 
the scale of unification of the GUT-normalised gauge couplings 

gi(M GUT )=g 2 (M GUT ). (1.4) 

We have applied the constraints at the gauge unification scale as an approximation. In 
principle, one should apply the constraints at the fundamental string scale, which may or 
may not coincide with the gauge unification scale. In the case that the string scale is not 
too many orders of magnitude greater than the gauge unification scale, our approximation 
should hold quite well. Later we will estimate this effect (or other possible corrections) by 
perturbing the fluxed boundary conditions imposed at Mqut- The strong gauge coupling 
is set from low energy data, and so ^(Mgc/t) / 9i,2{Mgut) is assumed to be shifted to the 
unified value by GUT scale threshold effects. We have checked that the required correction 
is at the per cent level only. 

We will use the following default inputs for Standard Model quantities, except where 
explicitly stated: the pole mass of the top quark filfl mj = 178 GeV, the running Standard 
Model value of the bottom quark mass [jl6[ m^m^)^ 1 ^ = 4.25 GeV, the strong coupling 
constant g§] a s {M z ) = 0.1187, «- 1 (M z ) M = 127.918 ||| and M z = 91.1187 GeV g 



2. Electroweak Symmetry Breaking 

The standard MSSM tree-level minimisation of the Higgs potential leads to the equations 



P tan 2 /3-l 2 z y ' 

fxB = - sin 2(3 (m 2 Hi + m\ 2 + 2fi 2 ). (2.2) 

As is well known, the minimisation conditions of the loop-corrected Higgs potential can 
also be cast in this form, with appropriate interpretations of the parameters (in particular, 
shifting fn 2 H . through tadpole terms). The minimisation conditions are imposed, as usual, 
at the scale Msusy = yJ m l Y m i 2 • ^ ne default version of S0FTSUSY1 . 9 fixes tan (3 and sgn(/j,) 
as input parameters and allows us to extract \^{Msusy)\ and B(Msusy)- In this approach, 



-m 



2 H2 tan 2 [3 + m 2 Hl 1 



2 This holds if fi is not specified. In an even more constrained scenario considered in p^ , in which both 
fi and B are predicted, dilaton dominance is not compatible with EWSB. 
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which we will follow in this section only, the value of B(Mgut) (i-e. re in eq. ( |1.3[) ) is 
computed rather than imposed. This will allow us to examine the consistency of electroweak 
symmetry breaking with the flux-induced soft terms. More specifically, S0FTSUSY1 . 9 allows 
us to interpolate the parameters of the MSSM below the weak scale and the GUT scale by 
integrating 2-loop MSSM renormalisation group equations (RGEs). It solves these RGEs 
while simultaneously imposing the boundary conditions eqs. Ql-l|), fll-4 ) an d adjusting the 
Yukawa and gauge couplings so that they agree with the data. For more details, refer to 
the SOFTSUSY manual @ • 




Figure 1: Electroweak symmetry breaking constraints in the fluxed MSSM for (a) \i > and 
(b) \i < and different values of M and to*. The horizontal line is the naive prediction of the 
model. The curves end when there is no viable model (see text). In (a), points are plotted where 
models exist which break electroweak symmetry correctly, but which possess a charged lightest 
supersymmetric particle (LSP). 

Fig. H shows the prediction for B(Mgut) as a function of tan (3, for \x > and \i < 0. 



The prediction is written in terms of re = B(Mgut)/{—2M). As eqs. (|1.2|) , ( [Oj ) state, the 
model predicts re = 1, which is indicated by a dotted horizontal line. We pick two different 
values of the mass scale M and three different values of mt and show the prediction for 
in each case. In Fig. |l|a we see that the correct value of rs for \i > is impossible 
to achieve for mt = 178 ± 4 GeV within 2a of its central value and for M = 300 GeV 
or M = 1200 GeV. We have checked numerically that this result holds for all M > 100 
GeV. The curves are truncated on the right hand side of the plot when < M^o , which 
is ruled out cosmologically for a stable neutralino lightest supersymmetric particle (LSP). 
For completeness, we also show extra points where such a condition is violated. The plot 
shows that significant departures from the r# = 1 prediction are required for the model to 
become viable (for example, tb = 0.85 can work for mt = 186 GeV). Fig. |l|b, on the other 
hand, shows that EWSB is viable for fx < and that there are two possible solutions for 
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tan (5, one around 3-5 and one at higher tan j3 > 25. This result holds for other values of M. 
The curves are truncated on the right hand side when m\ < 0, indicating an unsuccessful 
electroweak minimum. In summary, the plots in Fig. |] tell us that the model selects a 
particular sign of fi and two specific ranges of tan (3. 

Some features in Fig. |l| can be qualitatively understood in terms of the minimisation 
conditions and RGE effects. For instance, we can write B(Mgut) as a sum of two parts, 
B(Mgut) = B(Msusy) + AB, where B(M S usy) is determined by eqs. ( fOP , (fl2|), and 
AB encodes the RGE evolution between Msusy and Mgut, which reads 

Y rMGUT 

AB~—r (3h$A t + 3hlA b + h 2 T A T + g >2 M 1 + 3g 2 M 2 ) dlnQ (2.3) 

OTT JMsuSY 

at one-loop order. The AB contribution to B(Mgut) is always significant, and increases 
at large tan/3 because then the terms proportional to h 2 ,h 2 increase. On the other hand, 
eq. ( |2.2[ ) tells us that B (Msusy) is sizable at small tan/3 only, and that it flips sign if we 
switch from fj, > to fj, < 0. This explains why Fig. |l|a and Fig. [l]b are so different at 
small tan/3. At large tan/3 the plots are more similar to each other because B(Msusy) is 
suppressed, so the main contribution to B(Mgut) is AB, which is less sensitive to sgn(/j,). 
The dependence of AB on sgn(u) is mainly induced by hb,h T . These couplings are not 
exactly the same for /i > and /i < 0, because their relation with the fermion masses 
nit,, m T is affected by threshold corrections proportional to /xtan (3. As a result, Fig. |l|a and 
Fig. ||b exhibit some differences also at large tan /3. 

Also the fact that the curves in Fig. |l|a end before the corresponding curves in Fig. [l|b 
is mainly a consequence of the different values of hb,h T for positive or negative \x. As we 
have mentioned above, the curves end because either the lighter stau mass or the pseudo- 
scalar Higgs mass becomes too small (eventually tachyonic). The behaviour of such masses 
is shown in Fig. ||, for u, either positive or negative. Since the lightest stau is mainly f^j, 
the behaviour of ml reflects that of ml , as we have checked. The latter mass is M 2 
at Mgut- At smaller scales, m? fl receives negative RGE corrections proportional to h 2 , 
which are sizable at large tan /3 and eventually drive m 2 R to negative values. Stau mixing 
(proportional to tan/3) also decreases m^, but this is sub-leading to the RGE effect for 
the allowed regions of parameter space. As regards m 2 A , it is useful to recall its tree-level 
expression at large tan /3, which is m 2 A ~ rn/jj — m^ 2 — Mf . In fact, we have checked that 
the behaviour of m 2 A in Fig. ^ reflects that of m 2 Hi — m^ 2 . The latter quantity is driven to 
small values at large tan/3 because of RGE corrections proportional to h 2 or h 2 . 

The test of proper EWSB, which we have discussed above, is a crucial one for the 
viability of the model. Indeed, it is necessary for the existence of a phenomenologically 
realistic local minimum. Although this property is sufficient for the phenomenological 
discussions in the next sections, for completeness we should add that deeper minima seem 
to exist, with large Higgs vacuum expectation values (VEVs) and broken colour and/or 
charge. In particular, we have examined some field directions studied in [17], [18], which 
involve the Higgs doublet H2, the selectron doublet L e and a left-right pair of sbottom 
(or stau) fields. For sufficiently large field values, [f^l — \L e \ ^ M/hj (with hj = ht or 
h T ), the leading term of the potential along such a direction is V ~ ( m ^j 2 + m \ e )\H2\ 2 , 



- 5 - 



(a) 



m t =178 GeV, M 1/2 =300 GeV, u>0 
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Figure 2: Mass variation of the lightest stau, the pseudo-scalar Higgs and the lightest neutralino 
with tan/3 for M = 300 GeV, m t = 178 GeV and (a) [i > 0, (b) /x < 0. 



where the soft masses are evaluated at a renormalisation scale Q ~ | H2 \ ■ That sum of soft 
masses is positive at Mqut, then it decreases and eventually becomes negative, because 
of the behaviour of m^. This sign flip occurs at quite a large scale Q* ~ 10 9 GeV 3> M, 
hence the potential develops a very deep negative minimum in which \H2\ ~ Q* and 
\V\ ~ M 2 Ql 3> M 4 . The tunnelling rate between the realistic minimum and such a 
minimum is exponentially suppressed by a factor e~ s , where S is proportional to (or 
multiplied by a large numerical coefficient (see e.g. ]l9|). We will assume that this 
suppression is sufficient such that the probability of the universe to have tunnelled into 
the wrong minimum is small. In any case, a full computation of the tunnelling rate goes 
well beyond the scope of this paper. Note that the situation concerning charge/colour- 
breaking minima is similar to that in the case of another popular string-inspired scheme, 
dilaton dominated SUSY-breaking. Also in this case lower minima other than the standard 
EWSB minimum appear [13] and the assumption that we live in a metastable local vacuum 
provides a natural way out. 



In what follows, we use eq. ( |1.2[ ) (or eq. (|Q|)) as a boundary condition on B(Mgut)- 
We then predict tan (3 using two different numerical methods. If one imposes the tan (3 
prediction coming from eqs. ( |2.1| ), (|2,2j), it turns out that one always obtains the lower tan f3 
solution of EWSB. We therefore use this strategy when studying the low tan (3 solution. For 
the high tan/3 solution however, we revert to the default S0FTSUSY calculation, determining 
the correct value of input tan (3 by the method of bisection, using the fact that the curves 
in Fig. |l]b are monotonic near the high tan j3 solution. 
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3. Phenomenology 



We now detail the phenomenological constraints that we will place upon the model. We 
impose constraints coming from the decay b — > sj by calculating its value with the aid of 



micrOMEGAsl.3.1 |20|, |2l| linked to micrDMEGAS via the SUSY Les Houches Accord |2g]. The 
experimental value for the branching ratio of the process b — > is (3.52 ±0.30) x 10~ 4 p3| . 
Including theoretical errors [p4|| (0.30 x 10~ 4 ) coming from its prediction by adding the two 
uncertainties in quadrature, we impose 

2.3 x 10~ 4 < BR{b -> S7) < 4.8 x 10~ 4 (3.1) 

at the 3a level upon our prediction of BR(b — ► 57). micrOMEGAsl .3 . 1 also calculates the 
SUSY contribution to the anomalous magnetic moment of the muon, 5a^. The experimental 
measurement of a M = (g M — 2) /2 gives a very precise result, a£ v = (11659208±6) x 10~ 10 p5[ . 
It is difficult to predict the Standard Model value a^ M reliably at this level of accuracy. 
The estimates vary from being 0.7-3.2cr lower than the experimental number (see e.g. pfj). 



As a guide, we will take for a SM the value (11659189 ± 6) x 10~ 10 pi, which is 2.3a lower 



than a^ v . This would imply that the new physics contribution to is subject to the 3a 
constraint 

-6 x 10~ 10 < 5a^ < 44 x 10~ 10 . (3.2) 

Many authors (see for example f28[| ) have used the MSSM to explain the discrepancy 
between the experimental determination and the SM prediction of a^. It usually happens 
in the MSSM that 5a^ has the same sign as [i. This is true also in the special version of the 
MSSM we are discussing. In particular, since EWSB prefers negative [a, Sa^ is predicted 
to be negative, so the discrepancy is not alleviated in our scenario. Eqs. ( |3.1| ) and ( fT^ ) 
will prove to be strong constraints upon the fluxed \i < MSSM boundary conditions. 
The LEP2 collaborations have put stringent limits upon the lightest Higgs boson in the 



MSSM [29]. In the decoupling limit of Ma 3> Mz (applicable to the results discussed here 
where we apply limits to the lightest CP even Higgs boson), the 3a limit upon the mass is 
rrih > 114 GeV. The error upon the theoretical prediction is estimated to be typically ±3 



GeV for non-extreme MSSM parameters [3C] and so the bound 



m h > 111 GeV (3.3) 

is imposed upon the S0FTSUSY1 . 9 prediction of uih- 

We can also test the hypothesis that the lightest neutralino is a stable particle by 
examining its dark matter properties. We will compare the prediction of relic dark mat- 
ter density from thermal production in the early universe by micrOMEGAs with the 3 a 



WMAP Ipll , 32] constraint upon the relic density: 

0.084 < n DM h 2 < 0.138. (3.4) 

It is important to realise that the prediction of micrOMEGAs and the constraint from WMAP 
are made in the context of a 'standard' cosmological model (A— CDM). In the rest of the 
present paper, we will assume that the A-CDM model describes cosmology well. Then, any 
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prediction lower than the range in eq. ( |3.4|) would require additional forms of dark matter 
or non-thermal production and any prediction greater than that range would require the 
neutralino to be unstable (either if the neutralino were not the LSP or if R-parity were 
violated, for example) and for some other particle to constitute the dark matter. 

3.1 Low tan/? branch of \l < 

Fig. ||a shows the constraints coming from b — ► S7 and 5a^ for the low tan/3 branch of 
solutions to the fluxed SUSY breaking boundary conditions as a function of M, the free 
SUSY breaking mass scale. We see from Fig. |||a that the bounds coming from the 




Figure 3: Constraints on the low tan j3 branch coming from (a) BK(b — > sj) and the anomalous 
magnetic moment of the muon, (b) lightest CP-even Higgs mass and dark matter. The horizontal 
lines display (a) the lower bound on <5a M and the upper bound on BR(fe — » 57), (b) the lower bound 
on rrih coming from LEP2 constraints and the allowed region coming from the assumption that the 
neutralino makes up the entire dark matter relic density of the universe. 

anomalous magnetic moment of the muon and b — > S7 have a similar effect and imply 
that M > 270 GeV. From Fig. ||b, we observe that M > 550 GeV is required by the 
constraints upon the lightest Higgs mass m^. This last bound is notoriously sensitive to 



the value of mt (see e.g. [ 30 1 ) and will change significantly if we depart from the default 
central value of m t = 178 GeV. We will investigate this effect below, but we advance that 
the bound M > 550 is substantially relaxed and values M > 300 GeV are allowed for mt 
values 2 a above the central value mt = 178 GeV. The prediction of Qh 2 shows that the 
neutralino is not compatible with being stable dark matter in the low tan (3 branch. It is 
only compatible with the WMAP constraint (shown by the region between two horizontal 
lines on the figure) for M ~ 100 GeV, where the model is already excluded by m^, b — > S7 
and 5a M . 

The sparticle spectrum of the low tan/3 branch is shown in terms of M^o in Fig. |||. 
The lightest neutralino Xi is mainly bino, hence M o ~ 0.4M. From the figure, we see that 



-8- 




500 



(b) 



10 



7 

6 [ 

5 

4 

3 

2 



aW/h 4 



M(x?)/GeV 



100 200 300 400 
M(x?)/GeV 



500 



Figure 4: MSSM spectrum in the low tan/3 branch. The horizontal range corresponds to 100 
GeV< M < 1000 GeV. 

no other sparticles are particularly close in mass to the LSP, the nearest being the second- 
lightest neutralino and the lightest chargino. These two particles are quasi-degenerate (they 
are mainly winos) and are plotted as one curve on the plot, since individual curves could 
not be distinguished by eye. Also the second lightest chargino and the third and fourth 
lightest neutralinos are quasi-degenerate (they are mainly Higgsinos). In the slepton sector, 
the masses shown correspond to ex, £r, u e , but the same (or almost the same) results hold 
for the other generations. As regards lst/2nd generation squarks, the left-handed ones 
(not shown) are almost degenerate with the gluino, whereas the right-handed ones (qn) are 
slightly lighter. In the third generation sector we only show the stop squarks. The heavier 
sbottom, which is mainly 6r, is close to qn. The lighter sbottom, which is mainly 6x, is 
close to the heavier stop (tz), except for small M^o. The mass ratios do not vary much 
with increasing M except for the heavier stop. If the absolute masses were plotted, they 
would increase linearly with M (approximately). 

3.2 High tan (3 branch of \i < 

We now turn to the other branch of higher tan (3 that is valid for fj, < 0. Fig. |5|a shows 
that, contrary to the low tan (3 solution, Sa^ and b — * 57 provide much more stringent 
constraints: M > 750 and M > 1100 GeV respectively. This is expected, of course, 
because the SUSY contributions to both 5a^ and the b — > S7 amplitude grow with tan (3, 
hence a larger M is needed to provide the necessary suppression. On the other hand, Fig. |5|b 
shows that provides a much less stringent constraint of M > 140 GeV compared to 
the low tan/3 case. The main reason is that the tree-level contribution to m^, which is 
~ Mz \ cos 2/3\, increases for increasing tan/3. The neutralino is compatible with being a 
stable dark matter candidate for M ~ 200 GeV (which is already ruled out by the b — > S7 
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Figure 5: Constraints on the high tan j3, fi < branch, (a) constraints coming from the anomalous 
magnetic moment (with the horizontal line showing a lower bound) and from BR(b — > sj) (with 
the associated horizontal line showing an upper bound) . (b) constraints from the dark matter relic 
density and the lightest CP-even Higgs mass. Horizontal lines show the lower bound on mh, from 
the LEP2 direct searches and the WMAP constraints upon ilh 2 . 

and 5 dp constraints) or M ~ 1100 GeV, a remarkable point with a heavy MSSM spectrum 
that passes all constraints. 




Figure 6: MSSM spectrum in the high tan/3 branch. The horizontal range corresponds to 100 
GeV< M < 1300 GeV. 

Fig. |6] displays the sparticle spectrum in terms of M x o (which is about 0AM) for the 
high tan/? case. Some of the masses (neutralinos, charginos, gluino, 1st /2nd generation 
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sfermions) are not very different from those in the low tan (3 branch. Some other masses 
exhibit significant differences. In particular, each stau is lighter than the corresponding 
selectron 3 . Both sbottom squarks (not shown) have masses between t\ and t^. Finally, the 
heavy Higgs bosons are much lighter than gluinos and squarks, in contrast to the previous 
case. In particular, at higher values of M v o the value of approaches 2M v o. This 

_ Al Al 

enhances the x?Xi A ^ bb or rf annihilation channel through the pseudo-scalar Higgs 
pole, explaining the dip in £lh 2 displayed in Fig.5-b. 

As emphasised above, for M ~ 1100 GeV the lightest neutralino provides the appro- 
priate dark matter density and all other experimental constraints are fulfilled at the same 
time. The corresponding spectrum is quite heavy, with x? (the LSP) at around 500 GeV, 
fx around 600 GeV, X2 an d xf around 900 GeV, the heavy Higgs bosons around 1 TeV, 
and all other sparticles heavier than 1 TeV. This scenario could pose some problems on 
the experimental side, but SUSY should still be detectable at the LHC with sufficient lu- 
minosity. On the theory side, a possible draw-back is that such high M seems to require 
significant fine-tuning in EWSB, since Af| <C M 2 . For instance, proper EWSB requires 
|m|| <C M 2 at low energy, where m| = f^j 2 + A* 2 is the coefficient of l^l 2 - Hence an accu- 
rate cancellation is needed between the low-energy values of mn 2 H2 and fi 2 , which are a priori 
unrelated parameters with values 0(M 2 ). On the other hand it may be that some funda- 
mental model predicts a relationship between fi and the soft SUSY breaking parameters at 
the unification scale which leads to the appropriate cancellation at low energy. In such a 
scenario the fine-tuning problem would be alleviated, or even removed. In this respect, it is 
interesting to note that the values of /i(Mgut) required by proper EWSB are numerically 
close to — 2M, in our case. This is shown in Fig. |?| for both tan/3 branches. In particular, 
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Figure 7: The ratio h(Mgut)/M as required by EWSB, for the high and low tan/3 branches. In 
each case, the band is obtained by varying m t in its 2a range. The horizontal line corresponds to 
fi(M G uT)/M = -2. 



3 Sneutrino masses, which are not shown, are close to those of the corresponding charged sleptons (y e ~ e_L 
and v T ~ T2), except at small M x o. 
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we can see that a hypothetical relation (j,(Mgut) = — 2M is consistent with EWSB in the 
high tan (3 branch. It is also interesting that this happens for large M only (M > 900 
GeV), which is compatible with the experimental constraints discussed above. So, if such 
a relationship was found to hold due to some underlying physics at the unification scale 4 , 
one would have a model in which a little hierarchy Mz/M ~ 0.1 could emerge without 
requiring a special fine-tuning among mass parameters 5 . In such a model, in other words, 
the simple relationships among MSSM mass parameters at the unification scale, combined 
with the RG evolution of such parameters, would lead to cancellations in the low-energy 
Higgs potential, so that the smallness of Mz/M could be (at least in part) justified. 

3.3 fi > Case 

We turn now to the \i > set of solutions. Since we have found no solutions for the strict 
model prediction of r# = 1, we are forced to consider perturbations of the fluxed boundary 
conditions. Fig. |l]a demonstrates that higher rrit leads to values of rs closer to one, and 
that values of tb = 0.85 are attainable if rat = 186 GeV. We have taken these two values 
in order to study the scenario, but for purposes of brevity neglect to include the plots in 
this paper, choosing to summarise the results in text for M = 100 — 900 GeV. tan/3 is 
in the range 26-36, monotonically increasing with M. The LEP2 bound on m/j implies a 
mild constraint only (M > 150 GeV). The b — > 57 bound implies a stronger constraint 
(M > 500 GeV), due to the high value for tan/3. This bound, however, is weaker than 
the bound that holds for fj, < 0, because for \x > some cancellations take place in the 
b — > S7 amplitude. The advantage of [i > is even more evident in the case of Sa^, which 
is positive, so that the predicted a M more easily agrees with the measured value. In fact, 
the constraint from 8a ^ is not very restrictive (M > 230 GeV), despite the high values for 
tan p. 

The model's spectrum looks remarkably similar to Fig. |6[ The main differences are 
that the heavy Higgs bosons have larger masses, and the lightest stau is closer in mass to 
the neutralino. The latter feature results in two small regions in M where dark matter 
is compatible with the WMAP constraint. For the specified inputs, such regions have 
M = 160 - 190 GeV and M = 640 - 650 GeV. Between these two regions, Mi 2 is higher 
than the WMAP constraint, otherwise it is lower. 

4. Stability of Results 

We will now examine how stable the phenomenology is with respect to variations in exper- 
imental inputs, perturbations of the fluxed boundary conditions and theoretical uncertain- 
ties in the mass spectrum prediction. We will consider both tan/3 branches of the fj, < 

4 In fact it turns out that in simple toroidal orientifold Type IIB compactifications fj, and M correspond 
to different types of 3-form fluxes and in some simple cases they are related due to constraints coming 
from certain tadpole cancellations (see e.g.|3^, |34|). So it is not inconceivable that further constraints like 
\fl\ = 2M could appear in some specific situation. 

5 Only dimensionless parameters, i.e. gauge and Yukawa couplings, would need to be adjusted. This 
tuning, however, may be regarded as a separate problem. 
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case. In particular, an important question about the high tan/3 branch is whether it is 
possible to obtain a solution to all the constraints with a lighter spectrum than M ~ 1100 
GeV, taking into account various relevant uncertainties. 

In order to investigate the stability of our predictions, we will vary one input or bound- 
ary condition at a time, keeping all others at their defaults. These variations will be: 

• 2a variations in the input value of mt, i.e. 170 GeV < mt < 186 GeV. 

• 2cj variations in the input value of a s (Mz), i.e. 0.115 < a s (Mz) < 0.123. 

• Variations of MgusY by a factor of two in each direction (setting the minimum value 
to be not smaller than Mz), which should give an estimate of higher-order theoretical 
uncertainties (see e.g. ref. |30|]). 

• 10% variations in the boundary condition for M1/2 (0.9 < M1/2/M < 1.1), without 
changing the boundary conditions for uiq, Aq,B. 

• 10% variations in the boundary condition for B (that is, 0.9 < re < 1.1 in eq. (|1.3| )), 
without changing the boundary conditions for M 1 / 2 , mo > Aq ■ 

In practice, for a given value of M, we scan over 20 values of the parameter being varied, 
then plot the maximum and minimum values for the constraints obtained with that scan. 
Fig. |H shows the effect of the above variations on some of the predictions in the high tan /3 
branch (a,b,c) and in the low tan (3 branch (d). In each panel, we only show variations that 
have non-negligible effects. 

We first consider the high tan /3 branch. Fig. |8]a shows uncertainties in the relic density 
prediction. The uncertainty under variations of Msusy is small. Changing either the M1/2 
boundary condition by ±10% or the input value of a s (Mz) makes a larger difference, but 
does not qualitatively change the region compatible with WMAP, shown by the region 
between the horizontal lines. However, we see that varying either mt or the high-scale 
boundary condition on B makes any value 200 < M < 1200 GeV compatible with the 
WMAP constraint. As explained previously, above M > 200 GeV, low values of Qh 2 
are caused when two Xi s annihilate via almost-resonant s channel pseudo-scalar A Higgs 
bosons, i.e. for uia — 2M v o. Since m& is very sensitive to variations in mt or in the B 
boundary condition, the value of M at which the condition ~ 2M is satisfied can 

A 1 

shift very much under such variations. Hence the constraint on M from the relic density is 
not so stringent. However, Fig. |8|b shows that b — > 57 requires M to be anyhow larger than 
about 1 TeV, even if we allow for variations in mt or in the B boundary condition. Thus, it 
appears that we are stuck with a heavy spectrum (and the associated large fine-tuning and 
difficulties for detection in colliders) for the /i < high tan j3 branch. For completeness, 
we also show the behaviour of in Fig. [8|c, including variations. The constraint on M 
remains quite weak, but the absolute prediction on m^ has a strong dependence on mt, as 
usual. 

We turn now to the discussion of the low tan/3 branch. For default parameters, this 
branch is not compatible with the stable neutralino LSP assumption since the relic den- 
sity predicted is much higher than the WMAP constraint, as shown previously in Fig. |3|. 
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Figure 8: Uncertainties in predictions along the \i < high tan/3 branch in (a) the dark matter 
relic density fi/i 2 , (b) BR(b — > S7) (c) m^. Panel (d) shows in the fj, < low tan/3 branch. 
In (a), the region between the two horizontal lines is favoured by the 3cr WMAP constraint for a 
stable neutralino LSP making up all of the dark matter in the universe. 



Contrary to the case in the high tan/3 branch, this situation does not change with any 
variations, because none of them brings the model appreciably closer to a (co-) annihilation 
region. If we assume that this problem can be solved by some means (for instance by as- 
suming that the lightest neutralino is not stable on cosmological time scales), the strongest 
bound on M in the default calculation (Af > 550 GeV) comes from the LEP2 Higgs mass 
constraint, see Fig. [||b. Fig. ||d shows the behaviour of m,h under variations. We notice that 
rrih has the usual strong dependence on mt, as well as a significant dependence on the B 
boundary condition. Indeed, changing the latter parameter produces a significant relative 
change in tan/?, which in turn affects the tree-level Higgs mass (in this low tan/3 branch). 
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Fig. |8|i shows that varying either rat or the B boundary condition results in a lower bound 
M > 300 GeV. The variations induced by changes in Msusy or M]/ 2 have smaller effects 
upon the bound. The bounds coming from the anomalous magnetic moment of the muon 
and from b — * S7 do show smaller variations: the bounds on M can vary up to ~ ±10%. 



5. Conclusions 



We have considered the phenomenology of a set of SUSY-breaking MSSM boundary con- 
ditions motivated by flux-induced SUSY breaking in Type IIB orientifold models. These 
boundary conditions may be interpreted as the Type IIB version of modulus-dominance 
SUSY-breaking considered in the past for heterotic models. The model is extremely con- 
strained and, after imposing radiative EWSB, depends on a single parameter, the overall 
SUSY-breaking scale M. This is to be contrasted with other string inspired schemes like 
heterotic dilaton domination in which, even after imposing EWSB conditions, there are 
still two free parameters. We find it remarkable that consistent EWSB may be attained at 
all for our very restricted set of soft terms. 

We have shown that electroweak symmetry breaking is only compatible with one sign 
of [i < 0, for which there are two solutions: one with low tan (3 ~ 3 — 5 and one with high 
tan p~ 25 -40: 



The low tan (3 solution is consistent with having a SUSY breaking parameter M > 300 
GeV depending on the precise value of mt- This may allow for a relatively light SUSY 
spectrum. On the other hand this low tan (3 solution is not compatible with the light- 
est neutralino constituting the dark matter of the universe, since the predicted relic 
density is higher than that implied by WMAP and standard cosmology. Thus one 
would need to have an unstable neutralino, which could be due to the existence of an 
extra lighter SUSY particle or to the existence of R-parity violating couplings. In the 
latter case the collider phenomenology would very much depend on the neutralino 
lifetime. If it is much longer than O(10 _6 )s the lightest neutralino will leave detectors 
intact; collider signatures will therefore mimic the usual missing energy R-parity con- 



serving MSSM signatures and LHC discovery and measurements [35] should be easily 
possible. If the lifetime of the neutralino is shorter than 10 _6 s then the associated 
collider phenomenology will depend largely on its decay products. 

The high tan (3 branch of solutions can pass all constraints as well as provide the dark 
matter density compatible with WMAP constraints. The spectrum compatible with 
these constraints is heavy: with the LSP (x?) at around 500 GeV and all sparticles 
other than X2>Xi~ an d f\ heavier than 1 TeV. This would make sparticle measure- 
ments at the LHC practically impossible. With 100 fb _1 integrated luminosity, SUSY 
should still be detectable at the LHC through the inclusive i?™ lss measurement p6| , 
but more luminosity would be required for other channels involving leptons. 

Due to the heavy spectrum in this high tan (3 case a certain amount of fine-tuning of 
the EWSB conditions would seem to be required. On the other hand we have argued 
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that if the underlying theory predicted an approximate relationship fi = —2M at the 
unification scale, such fine-tuning would be substantially alleviated or even removed. 
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